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We suggest that white dwarf (WD) pulsars can compete with neutron star (NS) pulsars for produc- 
ing the excesses of cosmic ray electrons and positrons (e ) observed by the PAMELA, ATIC/PPB- 
BETS, Fermi and H.E.S.S experiments. A merger of two WDs leads to a rapidly spinning WD with 
a rotational energy (~ 10 50 erg) comparable to the NS case. The birth rate (~ lCP 2 -10 -3 /yr/galaxy) 
is also similar, providing the right energy budget for the cosmic ray e . Applying the NS theory, 
we suggest that the WD pulsars can in principle produce e up to ~ 10 TeV. In contrast to the NS 
\ model, the adiabatic and radiative energy losses of e are negligible since their injection continues 

after the expansion of the pulsar wind nebula, and hence it is enough that a fraction ~ 1% of WDs 
are magnetized (~ 10 7 -10 9 G) as observed. The long activity also increases the number of nearby 
sources (~ 100), which reduces the Poisson fluctuation in the flux. The WD pulsars could dominate 
the quickly cooling e ± above TeV energy as a second spectral bump or even surpass the NS pulsars 
in the observing energy range ~ lOGeV-lTeV, providing a background for the dark matter signals 
and a nice target for the future AMS-02, CALET and CTA experiment. 
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O ' I- INTRODUCTION 

h. 

Recently, the observational windows to the electron and positron (e*) cosmic rays are rapidly expanding the 

energy frontier, revealing new aspects of our Universe. The PAMELA satellite Q shows that the cosmic ray positron 

fraction (the ratio of positrons to electrons plus positrons) rises in the energy range of 10 to 100 GeV, contrary to 

the theoretical prediction of secondary positrons produced by hadronic cosmic rays interacting with the interstellar 

medium (ISM) [llOj . Shortly thereafter, ATIC/PPB-BETS @, suggest an sharp excess of the e ± with a peak at 

600 GeV, and although not confirming the ATIC/PPB-BETS sharp peak spectrum \ Fermi 0-11] and H.E.S.S 

, also suggest an excess of the e ± total flux around 100 GeV - 1 TeV compared to theoretical predictions based on 

t— I ■ low energy cosmic ray e ± spectrum (33l . |34| . All these observations of the excesses probably connected with the 

— ' PAMELA positron excess, and most likely su gges t a new source, possibly the astrophysical accelerators fl0l - [30l or 
<— ' ^ i i , , -i -i . • iSti nr^i i 



dark matter annihilation (3ll - [33l l3oT-f66ll /decay [36l l39l liil. [52l. loll. [6 4 6 ^ 83j] . although there might remain alternatives 



such as the propagation effects [89l493| or proton contamination~ [94l496l | . These discoveries have excited the entire 
particle and astrophysics communities and prompted over 300 papers within a year. See @ for a recent review. 

The most fascinating possibility for the e ± excesses is the dark matter, such as weakly interacting massive particles 
(WIMPs) that only appear beyond the Standard Model. Dark matter is a stable particle that accounts most of the 
matter in the Universe but the nature is not known yet. Usually, the observed excesses are far larger than expected 
in the conventional dark matter annihilation scenarios. The annihilation cross section must be enhanced by two or 
three orders of magnitudes larger than that for dark matter to leave the desired thermal relic density. Astrophysical 
boosts from substructure are difficult to accommodate such large enhancements. A possible solution is that dark 
matter interacts with a light force carrier, enhancing the annihilation by the Sommerfeld effect, only at the pres ent 
time (not at freeze out) [13, S3, H3l ■ The otherpossibilities include the dark matter decay [H H M, El H [MBlWIl 
and the annihilation boosted by resonances [H, H3] • Because the PAMELA anti-proton observations show no excess 
[8^.[85|. any dark matter model should preferentially produce leptons rather than hadrons. The other multi- messenger 
constraints with radio, gamma-ray and neutrino observations are also getting tight but not completely excluding the 
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1 The difference between the ATIC /PPB-BETS and Fermi results is still under debate |94H . In this paper we call these features as a whole 
" excesses" . 
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dark matter models [3J, |3J, |M IM IM 1^ 161 MM, 123, 129 IM _ 

More conservative candidates are the astrophysical accelerators in our Galaxy, such as neutron star (NS) pulsars fiol — 
[13] , supernova remnants (SNRs) p~8M26j[ , microquasars [28[ , or possibly a gamma-ray burst [H, H3] • Under plausibl e 
assumptions, they can supply sufficient energy for e ± cosmic rays, as already known befo re the PAMELA era [97l4ll4j ]. 
Cosmic ray propagate via diffusion in our Galaxy deflected by magnetic fields [ll5j | . Since cannot propagate 
far away due to energy losses by the synchrotron and inverse Compton emission, the sources should be located nearby 
(< f kpc). This proximity of the source provides a chance to directly probe the as-yet-unknown cosmic particle 
acceleration jf 1 2j and investigate how the e ± cosmic rays escape from the source to the ISM [27]]. Unlike dark matter, 



the astrophysical models generally predict, if at all, a broad spectral peak due to the finite source duration 10, 29]. 
The hadronic models such as SNRs also predicts the antiproton excess above ~ 100 GeV [H, H2] (but see 26|), as 
first pointed out by Fujita et al. [3, as well as the excesse s of secondary nuclei such as the boron-to-carbon and 
titanium-to-iron ratio [23], [25[ . The arrival anisotropy [29|, [98|, Il04j is also useful to discriminate between dark matter 
and astrophysical origins. The exc i ting thing is t hat these signatures will be soon proved by the next generation 
experiments, such as AMS-02 jll6l Ill7| . CALET jll8l Ill9l| onboard the Experiment Module of the International 
Space Station, and CTA |l2dl | on the ground, in coming several years. 

With the forthcoming next breakthrough, it is important to lay down the theoretical foundation for the TeV 
windows. In particular, there could still be room for additional astrophysical signals since the e cosmic rays have only 
< 1% energy budget of the hadronic cosmic rays. Although the supernova (SN)-related sources such as N S pulsars 
and SNRs may be the most plausible sources of the TeV e ± , there should be only a few local sources 12l| , while 
from distant sources can not reach us due to the fast inverse Compton and synchrotron cooling [l(J 112 j. Hence a 
clean window is possibly open for the dark matter or other astrophysical signals. A part of this window may have 
been already implied by t he spectra l cutoff around ~ 1 TeV in the H.E.S.S. data. The future AMS-02 experiment will 
detect e ± up to ~ 1 TeV |l!6t Ill7j. while C ALET will observe electrons up to ~ 10 TeV with an energy resolution 
better than a few % (> 100 GeV) [ll8llll9| . Also CTA will be able to measure the cosmic ray electron spectrum up 
to ~ 15 TeV (l20j . 

In this paper, we propose yet another source - white dwarf (WD) pulsars - that could potentially dominate the 
> TeV window or even already have been detected as the excesses above the conventional models [HI, [H| . A 
WD pulsar is an analogue of the NS pulsar with the central compact object being a WD rather than a NS. A spinning 
magnetize d compac t object generates huge electric fields (potential differences) in the magnetosphere via unipolar 
induction [122| - |l24l |. and accelerates particles to produce pairs if certain conditions are met. Then, almost all the 
spindown energy is transferred to the outflows of relativistic e^- , resulting in the cosmic ray . 

In our model, a rapidly spinning WD is mainly formed by a merger of two ordinary WDs (or possibly by an 
accretion), since the observed WDs are usually slow rotators [125] , Such a merger scenario was prop osed to explain 
Type la supernovae (SNIa). However, it is not clear that such mergers lead to the SN explosions [128| |. It seems 
reasona ble that a bout half of mergers leave rapidly spinning WDs with the event rate of about one per cent ury i n ou r 
Galaxy [HI [l27|- The strong magnetic fields (> 10 6 G) are also expected as a fraction - 10% of WDs (l2 a ll3C j . 
Combining these facts, we will estimate that the WD pulsars can potentially provide the right amount of energy for 
the cosmic ray e ± (see Sec|TT|. We note that the WD mergers are also related to the low frequency gravitational wave 
background for LISA 

The WD pulsa rs have been theoretically adopted to interpret the observational features of the anoma lous X-ray 
pulsars |l32| - jl34j |. the close binary AE Aquarii [135| . and the transient radio sou rce GCRT J1745- 3009 137 1 . Our 



calculations for the e^- production are essentially similar to that of Usov |l33l . Il34| and Zhang & Gil [l37j |. However, 
this is the first time to apply the WD pulsars to the cosmic rays, as far as we know. We also discuss the adiabatic 
energy losses of e^- in the pulsar wind nebula, that are found to be negligible in contrast to the NS model. From 
the observational viewpoint, the WD pulsars have not been firmly e stabl ished, whereas there are several indications 
for their existence, such as the hard X-ray pulsation in AE Aquarii [l36j . The WD pulsars are likely still below the 
current level of detection because they are rare, ~ 10~ 4 of all WDs, and relatively dim. 

This paper is organized as follows. In Sec|Hl we show that the WD pulsars can produce and accelerate e^- up to the 
energy above TeV. At first we show that the energy budgets of WD pulsars are large enough to explain the PAMELA 
positron excess by order-of-magnitude estimates. Then we discuss, more closely, whether or not WD pulsars can 
produce and accelerate e up to the energy above TeV by considering the magnetospheres and pulsar wind nebulae. 
We also point out that there should be much more nearby active WD pulsars compared with NS pulsars since the 
lifetime of WD pulsars are much longer. In Sec lIIH we discuss the propagation of the from WD pulsars, and show 
the possible energy spectrum observed by the current and future observations in the WD pulsar dominant model 
and the WD and NS pulsar mixed model. As complements, we also give a short review of the current status of the 
observations of WD pulsar candidates. In Sec lIVl we summarize our paper and discuss open issues. 



3 



II. WHITE DWARF PULSARS 
A. Energy Budgets of White Dwarf Pulsars 

In this subsection, we show that WDs potentially have enough rotational energy for producing high energy e^- 
cosmic rays. 

NS pulsars, which are formed after the SN explosions, are one of the most promising candidates for the astrophysical 
sources of high energy positrons. For the PAMELA positron excess, each NS pulsar should provide mean energy 
~ 10 48 erg to positrons [HI [ll|, since the energy budgets of cosmic ray positrons is ~ 0.1% of that of cosmic ray 
protons, which is estimated as ~ 10 50 erg per each SN, and the positrons suffer from the radiative cooling during the 
propagation more than the protons. The intrinsic energy source is the rotational energy of a newborn NS, which is 
typically 

" > ~ m '° {i£k) («*=)' (toS-t) 2 ^, ID 

where I is the moment of inertia of the NS. Then, if all the NS pulsars are born with the above rotational energy 
and the ~ 1% energy is used for producing and accelerating , the NS pulsars can supply enough amounts of for 
explaining the PAMELA positron excess |10j . 

Let us show that double degenerate WD binary mergers can also supply enough amounts of rotational e nergy . Here 
we consider the mass 0.6M Q and radius R ~ 10 8 ' 7 cm for each WDs, which are typically observed ones |l4l| . Just 
after a merger of the binary, the rotational speed v ro t can be estimated as v to t ~ (GM/R) 1 / 2 ~ 10 8 cm/s, which 
corresponds to the mass shedding limit, and the angular frequency is about ft — v ro t/R ~ 0.1s -1 . Then, the rotation 
energy of the merged object is 

**™ - ~ ios ° Gab) g«=)' (o^y^ <2 > 

which is comparable to the NS pulsar case in Eq. ([1]) . The event rate ??wd of the double degenerate WD mergers in our 
Galaxy remains uncertain. Any theoretical estimate requires a knowledge of the initial mass function for binary stars, 
the distribution of their initial separation, and also the evoluti on of the system during periods of nonconservative 
mass transfer. There are still reasonable estimates in the range [l26l Il27l |. 

WD ~ 10~ 2 -10- 3 /yr/galaxy. (3) 

This is comparable to the typical birth rate of NS pulsars [l42l Il43| . Therefore, from the viewpoint of energy budget 
in Eqs. {TJ), and ©, the WDs are also good candidates for the high energy e ± sources as the NS pulsars, if the 
merged binaries can efficiently produce and accelerate e ± . 

The estimated merger rate is also similar to that of SNIa, which is one of the reason that the double degenerate 
WD mergers are possible candidates for SNIa. Since the typical WD mass is 0.6M Q , the merged objects do not exceed 
the Chandrasekhar limit lA M Pl even without any mass loss. Then, they leave fast rotating WDs, as suggested by 
some recent simulations [l44l |. and could become WD pulsars. In this paper, we assume that a fair fraction of double 
degenerate WD mergers result in the WD pulsars. 2 

The accretion scenario is another possibility for the fastly rotating WD formation. In the single degenerate binary, 
which consists of a WD and a main sequence star, there should be a mass transfer from the main sequence star to 
the WD as the binary separation becomes smaller and the Roche radius becomes larger than the radius of the main 
sequence star. In this stage, the angular momentum is also transferred to the WD, and the WD can spin up to around 
the mass shedding limit with the rotational energy as large as Eq.©. In Scc lIIIDl we refer to such a WD pulsar 
candidate, AE Aquarii. 

Since the birth rate is relatively uncertain in the accretion scenario, we just concentrate on the merger scenario in 
this paper. 



2 Since the highly magnetized WDs have higher mean mass ~ O.95M0 than the total average ~ 0.6Mq I129H . the fraction of mergers that 
leave spinning WDs could be lower than the average. 
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B. e Production and Acceleration 

In this subsection we discuss the possibility that WD pulsars emit high energy above TeV. In order to produce 
the TeV e^, a pulsar has to 

(i) produce pairs 

(ii) accelerate up to TeV. 

We show that WD pulsars can meet both of the conditions. From now on we set fiducial parameters of the WD pulsar's 
surface dipole magnetic field, angular frequency, and radius as B p = 10 8 G, SI = O.ls^ 1 and R = 10 8 ' 7 cm, respectively. 
For comparison, we set fiducial parameters of the NS pulsars as B p = 10 12 G, ft — 10 2 s _1 and R = 10 6 cm. 



1. pair production in magnetosphere 

Some of the observed WDs have strong mag netic fields of B ~ 10 7 " 9 G 129, 130]. For such WDs, if they are rapidly 
rotating as we discuss in the previous subsection, the electric field along the magnetic field are induced on the surface 
and the charged particles are coming out from the surface layer of the pulsars. Then we can expect that, as in the case 
of ordinary NS pulsars, the corotating magnetosphere are formed around th e WD s, in which the charge distribution 
of plasma should be the Goldreichi- Julian (GJ) density in a stationary case |l22| . 

(n x r) x B f2 B 10 5 / B p \ ( fl \ „ 

p0 = v - Jc *-^~-W\K^G){^) cm ' (4) 

where Z is the elementary charge of particles in the plasma. Here we assume that the large scale configuration of the 
magnetic field is dipole. Since the corotating speed of the magnetic field lines cannot exceed the speed of light, the 
magnetic field cannot be closed outside the light cylinder R\ c = c/fi. This fact lea ds to the open magnetic field lines 
in the polar region. The electric potential difference across this open field lines is [l22j [ 

_ b p &r* / Bp w a \ 2 / r \\ 

AKmax_ 2 C 2 10 ^108gJUib-V ll0 87 cmj Volt ' (5) 

which is the maximum value for the pulsars in principle. 

If the GJ density is completely realized in the magnetosphere, electric fields along the magnetic field lines is absent: 
E • B = . Since the charged particles are tied to the strong magnetic field, the acceleration of e ± cannot occur. That 
leads to the absence of high energy 7 ray emissions from the accelerated e^ 1 and successive pair production avalanches. 
However, there are two prospective scenarios of forming the region where the charge density is not equal to the GJ 
density, and hence are accelerated and produced in the NS pulsar magnetosphere, that is the polar cap jl23l Il45| 
and outer gap model |l24j . From now on, we assume that the magnetosphere structure of WD pulsars are similar to 
that of NS pulsars, and discuss the e 1 *" pair production especially in the polar cap region. 

In polar cap models, electric potential drops along the magnetic fields are formed in the polar region of the pulsars. 
There are some different types of polar cap models. First, the angle between the magnetic and rotational axis 
determ ines the sign of electric charge of the particles propagating along the open magnetic field lines in accordance to 
Eq.dU) [122( | . The GJ density in the polar cap region is positive when f2 • B < and negative when J~2 ■ B > 0. Second, 
polar cap models depend on whether or not steady charge currents flow out from the surface of the pole region. After 
the GJ density is realized, there are no electric forces working on the charged particles in the surface layer. Hence, 
whether or not the charged particles come out from the surface is determined by the competition between the binding 
energy of i ons o r electrons at the surface and thermal energy. In the original model proposed by Ruderman and 
Sutherland jl23| . they assume that the binding energy is bigger. Then due to the outflow along the open magnetic 
field, a gap where the charge density is almost is formed in the pole re gion. On the other hand, if the thermal 
energy is bigger, there exist a positive or negative space-charge-limited flow |l45| |. Even in this case, it is shown that, 
by virtue of the curvature of magnetic fields, the charge dens ity deviates from the GJ density and electric potential 
drops along the open magnetic field lines can be formed |l45j . Although a g eneral relativistic frame dragging effect 
also contributes to form electric potential drops in the polar cap regio n |l46| |. the effect can be neglected compared 



with the effect of magnetic field curvature in the case of the WDs 137]. 

In the polar cap region, where the GJ density is not realized, primary electrons or positrons are accelerated, and 
they em it cu rvature radiations, which interact with the magnetic fields and produce secondary e ± pairs, 7 + B — > 
e~ +e + jl23l |. The secondary e ± are also accelerated and emit curvature radiations that produce further e ± pairs (pair 
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creation avalanche). Inv erse Compton scatterings can also serve as a way to produce high energy e ± and successive 
pair creation avalanches [l47j . Due to the abundant charges supplied by the avalanche, the GJ density is realized at a 
finite distance from the surface and the polar cap formation stops. In the quasi-steady state, the size of the polar cap 
region can be approximated as h re I, where I is the mean free path of the pair creation process. To put it the other 
way around, only when the available size of the polar cap region /i max is larger than I, e pair creation avalanches 
can be forme d. Ch en & Ruderman first derived the condition for NS pulsars and succeeded in showing the NS pulsa r 
" death line" jl48| | . Harding & Muslimov also derived the NS death line under more general conditions jl70l Il7ll | . 
Here we follow Chen & Ruderman's approach and drive the e^- pair production avalanche condition in the case of 
WD pulsars. We discuss the validity of this simple treatment in Sec lIVI 

Going through any potential drop AV along the open magnetic field lines, are accelerated up to the Lorentz 
factor 

7 = ZTT2 . ( 6 ) 



where m c is the mass of electrons. The characteristic frequency of curvature radiation photons from the accelerated 



is 



7 3 -, (7) 



where r c is the curvature radius of the magnetic field lines. The mean free path of a photon of energy hui > 2m e c 2 
moving through a region of magnetic fields is jl49l | 

He h B a ( 4 \ , 

I = 4.4 -2L exp — ; X « 1), 

e z m e c B± \^X/ 



= 1}u B ± 

X ~ 2m e c 2 B q ( ' 

^ q — i,l c ^ /eh = 4.4 x 10 13 G and B± = B s sin9 with 8 is the angle between the direction of propagation for 
photon and the surface magnetic field lines of the pulsars. B s is the local magnetic field at the surface of the pulsar 
which is not necessarily coincident with the dipole field B p . At distance h above the pulsar surface, the sin 9 can be 
approximated to re h/r c , then 

B±_ re B s -. (9) 

r c 

We shall consider the situation I re h, however which could be realized when = O(10) without relying on precise 
value of the parameters characterizing NS or WD pulsars since small c hang es in x correspond to the exponentially 
large change in I. Here we take the critical value as x 1 = 15 following [148]. Substituting Eq.©, ©, (JSJ) to Eq.©, 
this condition is given by 



eAV\ 3 h h B s 1 



m e c 2 J 2m e cr c r c B„ 15 



(10) 



Eq. flTO]) corresponds to a general condition for e ± pair production avalanches in the polar cap region of pulsars. Then 
we have to specify h, AV, B s and r c . The thickness h and the potential drop Ay in the polar cap region depend 
on which p olar cap model we adopt. Here we consider the original polar cap model proposed by Ruderman and 
Sutherland |l23l j . In this case, the relation between h and AV is given by 

B s VLh 2 , . 

AV = — . 11 

2c v ; 

Then, since AV cannot exceed the maximum potential drop available in a pulsar magnetosphere, Ay max in Eq.©, h 
also cannot exceed the maximum thickness 

h ra «(^) 1/2 . (12) 
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B s and r c depend on the configuration of the surface magnetic held, which is very uncertain even in the case of the 
NS pulsars. Here we suppose curved magnetic helds in the polar cap region and set r c rs R and B s rs B p . In this 
case, the condition for pair production avalanche (Eq.([T0|)) is 

/ eAVmax \ 3 h femax B p > 1 

\ m c c 2 J 2m c cR R B q ~ 15' ^ ' 

which is equivalent to 

4 log B p - 6.5 log P + 9.5 log R > 96.7, (14) 

where the unit of B p , P = 2n/fl and R are [G], [sec] and [cm], respectively. By substituting R ~ 10 6 cm, which is the 
typical radius of NSs, Chen and Ruderman succeeded in explaining the NS pulsar death line [1481 ] . In the case of WD 
pulsars, substituting our hducial parameters B p ~ 10 8 G, P ~ 50s (CI ~ 0.1s) and i? ~ 10 8 7 cm, we find that the WD 
pulsars well satisfy Eq. ([T4l . and thus also the condition (i) in Scc lII Bl 

FigQ] shows the death lines of the WD and NS pulsar with the fiducial parameters of the WD pulsars. We also 
plot parameters of the observed WD pulsar candidates, AE Aquarii and EUVE J0317. As we discuss in Sec lIII DI the 
pulse emission like ordinary NS pulsars are observed for AE Aquarii, and not for EUVE J0317, which is consistent 
with the death line. 




FIG. 1: This figure shows the death lines of WD (solid line) and NS (dashed line) pulsars. The cross shape indicates the fiducial 
parameters of the WD pulsars, B p — 10 8 G and P — 50s. The observed data of rapidly rotating magentised WDs, AE Aquarii 
(filled circle) and EUVE J0317 855 (open square) are also plotted. The parameters and observational properties of these WDs 
are given in Sec lIII Dl 



2. e acceleration and cooling in pulsar wind nebula 

In the previous subsection we show that WD pulsars can produce pairs in the magnetospheres. In this subsection 
we discuss the acceleration and cooling of the in the pulsar wind nebulae. 

FigH] shows the schematic picture of a expected WD pulsar wind nebula. Once a WD pulsar is formed, the 
relativistic wind blasts off from the pulsar magnetosphere ~ i?i c . The supersonic wind becomes subsonic by passing 
the shock front at ~ Ri n , reaches the ISM and forms a contact discontinuity. Since the wind is continuously injected 
by the pulsar, the contact discontinuity keeps sweeping the interstellar matter, and then the outer shock front is 
formed at ~ i? ut- We emphasize that the SN shock front does not exist outside the shocked region unlike the NS 
pulsars since there suppose to be no SN explosion when the WD pulsar is formed. 

First we estimate the energy of available in the wind region R\ c < r < _R; n . In principle, can be accelerated 
to the energy that the equipartition is realized between the wind and magnetic held, eN — B 2 /8tt , that is 



B 2 

8ttV 



(15) 



where N is the number density of e . If the number flux is conserved in the wind region, 4nr 2 cN w const, N can be 
described as 

N = N lc (^\ , (16) 
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FIG. 2: This figure shows the schematic picture of the expected WD pulsar wind nebula. 



where Ni c is the number density at the light cylinder which can be estimated as 

Mc = ^M = f^V (17) 
e Zirce 

where p\ c and B\ c are the GJ density (Eq.Q) and magnetic field strength at the light cylinder, respectively and M. 
is the multiplicity of e^- in the magnetosphere. Inside the light cylinder r < R\ c = c/O, the magnetic field is almost 
pure dipole, 



For the fiducial parameters of WD pulsars, the radius of the light cylinder is R\ c ~ 3 x 10 11 cm and B\ c = 
B P (R/ (c/il)) 3 ~ 1G. Outside the light cylinder r > Ri c , if the energy flux of the magnetic field is also conserved 
B ■ r « const, then 

B = B lc — • (19) 
r 

Substituting B\ c = B p (QR/c) 3 and Ea flTTj) to Ea. (|T5|) . the typical energy of in the wind region can be described 
as 

where AV max is shown in Eq.([5]). The multiplicity of e in the pulsar magnetosphere and wind nebula have not been 
understood clearly even in the case of NS pulsars and there are several discussions |lOll . Il02j | . Although details of 
the multiplicity in the magnetosphere cannot be discussed at this stage 3 , TeV energy could come out of the wind 
region and the condition (ii) in Sec lIIBI can be fulfilled if M. is not large. 

Secondly we estimate the adiabatic and radiative cooling of e^- in the shocked region. To that end, we have to 
identify the radii of the inner and outer shock front i?i n and i? ut- The equation of motion for the outer shock front is 



dt 1 3 outP dt 



<j —Rl utP —^l \ = 4^R out 2 P sh , (21) 



3 In [l33j . Usov discussed the multiplicity in the magnetosphere for a X-ray pulsar IE 2259+586 based on the WD pulsar model by 
investigating the observed X-ray luminosity, in which M ~ 0.1. 
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where P s h is the pressure of the shocked region and p is the density of the ISM p ~ 10 24 g cm 3 . The energy 
conservation law at the outer shock front is 

d (4:71 ,3 1 d f 4tt ,1 

*{t^ V^^-^ItM' (22) 

Here L is the spin down luminosity of WD pulsars, 

BlQ 4 R e 

L = S—, (23) 

and we suppose that in the shocked region the particles are relativistic and its internal energy is 3P/2. Solving Eq. (|2"Tj) 
and (J22J) for R out (t), 



1/5 / r \ 1/5 

1547T J 
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16 



{ B p \ 2/5 { Q \ 4/5 / P \ 6/5 / t^ 3/5 l " 4 ' 



li0 8 G/ I O.ls" 1 / ll0 8 7 cm 



yr 



The outer shock finally decays when the pressure of the shocked region P s h becomes equal to that of the ISM p. At 
this stage the shocked region may be physically continuous to the ISM. Solving Eq. (|2~Tj) and (|2"2")l for P s h, 

- io (iwj (oi^j (to^j y dyn/cm ■ 

Besides assuming that the density of the ISM is p ~ 10~ 24 g cm~ 3 , that is the number density of hydrogen is n ~ 1cm -3 , 
the pressure can be estimated as 

13 ( T \ j„„ /„ m -2 



p = nfc B T ~ 10- 13 (^Y^J dyn/cm \ (26) 

where fee = 1.4 X 10 _16 erg K~ is the Boltzmann constant and T is the temperature of the ISM. From Eg. (1231) and 
(|26|) . the outer shock decays at about 

/ T \ 5/4 

* dec ~ 10 [wk) y1 "' (27) 
for the fiducial parameters of the WD pulsars. The lifetime of a pulsar r can be estimated as 

Prot 

From Eq.([2]) and (|23|) . for the fiducial parameters of the WD pulsars 



(28) 



Compared with Eq. ([2T|) and (|21)1) . we found that the outer shock decays at a very early stage of the lifetime of WD 
pulsars. 

For t < tdec, the momentum transfer by the wind balances the pressure of shocked region at the inner shock front, 

Lc , . 
2 = Ph. (30) 



4ttP; 

Then the radius of the inner shock front can be estimated as 



1/2 i/5 / T x 3/10 

(31) 



, 2/5 

]0 15 (-?-J cm, 



9 



for the fiducial parameters. For t > idee, there is no well-defined shocked region any more and the radius of the inner 
shock front is determined by the balance between the wind pressure and the pressure of the ISM p instead of P sn , and 
Ri n (t) become constant for t. For the fiducial parameters, 

Rin(t > t dcc ) ~ 10 17 cm. (32) 

In the case of NS pulsars, the adiabatic cooling due to the expansion of the shocked region is considerable as a cooling 
process in the pulsar wind nebula. However, in the case of WD pulsars, since the outer edge of the shocked region 
does not expand after t > idec, the adiabatic cooling shall give minor contributions to the cooling process of the high 
energy e . 

Now we discuss the e ± radiative cooling in the shocked region r > R m . In the region swept by the shock, the 
magnetic field may be highly fluctuated and the high energy e coming from the wind region are trapped because of 
the multiple scattering by the field, and lose the energy by the synchrotron radiation and inverse Compton scattering. 
Here we take the Bohm limit, where the fluctuation of the magnetic field SB is comparable to the coherent magnetic 
field strength B. In this limit, the diffusion coefficient -D s h can be approximated by 

Ah = (33) 

where r g = e/eB is the Larmor radius of the with energy e. The time scale t<af for the trapping in the shocked 
region is given by 

d 2 3eBd 2 

£dif = = ^ , (34) 

where d is the size of the shocked region. 

We consider the age t = t\vd > £dcc- For t > idec, we set the size of the shocked region as the forward shock front 
at t = idee, that is 

d « flout (t = W) ~ 10 19 cm, (35) 

for the fiducial parameters. As we have shown in Eq.(|32p. the radius of the inner shock front is about Ri n ~ 10 17 cm 
at t = t\vd- From Eq.(|19p. the strength of the magnetic field at the inner edge B ln can be estimated as 

Rh 



A.-3X10 "M lo^r) G ' (36) 



which is almost the same as that of the ISM. Then, substituting Eq.flSS]) and Eq.® to Eq.([3"5|) 4 , the time scale for 
the high energy e^- with energy e being trapped in the shocked region is 

-1 

,10TeV> 

The synchrotron energy loss of the e ± with energy e is described as 



<dif ~ 3 x wi (wkv) yr - (38) 



— - --a c/3 2 — ( C ) (39) 
dt 3 8ir \ m c c 2 J 

where ctt is the Thomson scattering cross section, and j3 — v c /c is the velocity in terms of the speed of light. Then 
from Eg. (1391) . the typical energy loss of the electron with energy e during the time scale tdif can be estimated as, 

This means that the high energy injected into the shocked region lose roughly 10% of the energy by the synchrotron 
radiation before diffusing out into the ISM. Although the inverse Compton scattering is also considerable process as a 
radiative cooling, it would be comparable to the synchrotron cooling. Then we can conclude that the radiative energy 
loss of in the pulsar wind nebula is not so large. 



4 In this case, the diffusion coefficient can be estimated as 



D sh - 10 : 



24 f 



; ) cm 2 /s. (37) 



V3GeV' 

This D s h is smaller than the diffusion coefficient in the ISM (see Eq.{5Sj) , which means that we consider the situation where the are 
highly trapped in the shocked region. 
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C. Differences between white dwarf and neutron star pulsars 

In this subsection, we discuss the differences between WD pulsars and NS pulsars as TeV e^- sources. 

Ordinary NS puls ars have been already discussed as a candidate for high energy sources for the PAMELA 
positron excess ([10| and the references listed in Secfl]). Compared with the NS pulsars, there are distinct features of 
the WD pulsars as high energy e ± sources. As we saw in the previous sections, the WD pulsars can provide the high 
energy and the intrinsic energy budgets are almost the same as that of the NS pulsars. However, the magnetic 
field and rotation speed of the WD pulsars are much smaller than that of the NS pulsars. As a result, the spin down 
luminosity (Eq.(|23p) of the WDs are much smaller than that of the NSs, 



wd ~ 10 41 (") ) erg/yr ~ 10"^ NS . . 41 ) 



6 

\Wg J Vo.is- 1 J V 108 - 7cm 

Then from Eq.(|28p. the lifetime of the WD pulsars are much longer than the NS pulsars 

twd ~ 10 9 yr~ 10 4 t ns . (42) 

Therefore, the number of the WD pulsars which are currently TeV sources are much larger than that of the NS 
pulsars. Since the high energy electrons above TeV cannot propagate more than ~ 1 kpc in our Galaxy, the number 
density of the WD pulsars which can be the TeV sources is 



a ■ ?7wd • t W d ^ 1q3 / WD \ f twd \ f Vg 



nWD = V G 10 "U-V-galaxy-J [W&) [To^ ) ^ " (43) 

where Vg is the volume of our Galaxy and ?7wd is the event rate of the double degenerate WD binary merger in our 
Galaxy, Eq. © . A parameter a is the fraction of the binary mergers which lead to the WD pulsars with the strong 
magnetic field B > 10 8 G. Eq. (|43t means that there may be enough WD pulsars which supply TeV near the Earth, 
although the parameter a has a large ambiguity at this stage. On the other hand, the number density of the TeV 
sources for the NS pulsars is 

n NS ~ O.lkpc" 3 - IO^oT^wd- (44) 

Eg. (1441) means that it is uncertain whether NS pulsars are sources above TeV energy or not. 

Another important difference is the environment of the pulsars, especially the strength of the magnetic field in the 
pulsar wind nebulae. The magnetic field is crucial for the cooling process since it determines how the high energy e^- 
produced at pulsars are trapped and lose their energy by synchrotron radiation in the pulsar wind nebulae. In the case 
of the WD pulsars, the strength of the magnetic field at the shocked region is, in most of their lifetime, comparable 
to that of the ISM. As we saw in the previous subsection, this may imply that most of the accelerated e^- directly 
escape into the ISM without cooling in the shocked region. On the other hand, in the NS pulsar wind nebulae, the 
situation is quite different. First the magnetic field are much stronger than the WD pulsars. Second there exists a 
SN shock front outside the pulsar wind nebula. These facts make the cooling process in the pulsar wind nebula more 
complicated, and the escape process of into the ISM is still uncertain. 

In the case of the NS pulsars, almost all the spin down luminosity is transformed to the kinetic energy of the 
e^- wind before the wind goes into the shocked region [l50j ]. The NSpulsars are consistent to be the source of the 
observed if the e^- lose ~ 99% of their energy in the shocked region 101. As we discussed in Sec|Hl the total energy 
budgets of the WD and NS pulsars are almost the same when almost all the double degenerate WD binaries merge 
to become the WD pulsars, that is when a = 1. Since the lose only ~ 10% of their energy in the WD pulsar wind 
nebulae fEq. (l4T)l) L the expected amount of from WD pulsars can exceed the current observation bound. Hence if 

a ~ 0.01, (45) 

we can expect that the PAMELA positron excess can be explained by the WD pulsars without any contribution 
of other sources. We should note that the fraction in Eq. pBI) seems consistent with the observed fraction of the 
magnetized WDs ~ 10%. We show the brief summary of the comparison between WD and NS pulsars in Table HI 



III. ENERGY SPECTRUM CALCULATION 

In this section, we calculate the energy spectrum observed at the solar system after the propagation in our 
Galaxy for the WD pulsar model. We solve the diffusion equation taking into account the Klein- Nishina (KN) effect. 
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pnprfxv npr papVi 

Evot [erg] 


1 1 1 in i n oqi t" v npr pspVi 

L LlllllllWol \j V l_/*_.± rl^cVL^ll 

L [erg/yr] 


lifetime 
T W 


event rate 
[1/yr/galaxy] 


number density 
n [1/kpc 3 ] 


pffi P1PT1PV 

[%] ' 


L x n x efficiency 
[erg/yr/kpc 3 ] 


WD pulsar 


~10 50 


~10 41 


~10 9 


~ a/100 


~ 10 3 a 


~ 90 


~ 10 44 a 


NS pulsar 


~10 50 


~10 45 


~ 10 5 


~ 1/100 


~ 0.1 


~ 1 


~10 42 



TABLE I: The comparison between WD and NS pulsars as e sources. 



A. electron distribution function from a single source 

Here we formulate cosmic ray propagation through our Galaxy according to [l07j . For simplicity we assume that 
the diffusion approximation is good (e.g., neglecting convection), the propagate in spherically symmetric way and 
diffuse homogeneously in our Galaxy 5 . Following these assumptions, the propagation equation can be described 
as follow. 

Here /(t, e,r) [m~ 3 • GeV~ x ] is the energy distribution function of e^. P(e) is the cooling function of the e^- which 
corresponds to the energy loss rate during the propagation. D(e) denotes the diffusion coefficient, which does not 
depend on the position r. Q(t, e, r) is the energy injection term. Considering 5-function injection at the time t = to, 
that is 

Q(t,e,r) = AN(e)d(r)6(t-t ), (47) 

we can get the analytical solution \lOT\ . For an arbitrary injection spectrum AN(e), the energy distribution can be 
described as 

/(r , t , e ) = ^f4^) e xp(-4). ,48, 
n 6 i z r A if "(e) V r dif / 

Here tt,o corresponds to the energy of which are cooled down to e during the time i — to, and is obtained by solving 
the integral equation 

r t0 de' , . 

t~to= ^rr. (49) 



P(e'Y 



The e propagate to the diffusion length defined by 



/ ftt.o D ( e >) \V2 

r dif (e,e t , )=2^ T(^ de ') " ( 50 ) 

Eg. (1481) is the distribution function for the S- functional (short term) injection source, i.e., the Green function of 
Eq. (l46l) . From now on, we set the observation is taking place at t = 0. 

Even for a continuous (long term) injection source, the distribution function can be calculated by integrating Eq.(|48j) 
for the active time of the source. The integration can be done numerically by transforming the integration from dto 
to de t ,o = P(e t fi)dt . (That is, we take e t ,o a s the time coordinate.) Substituting AiV(e tj o(e, to)) = Q(et,o(e,to))dto 
into Eq.(|48p and integrating over dto, the resulting distribution function reads 

J^ r > = 3/2 VI \ / 1 U eXP 1 ^2 j de t,0- (51 

7r J /^F(e) J e r d if(e, e t ,o) d \ r di f(e,e t )V 
Here e t - is the energy of e ± when they leave the source at the source birth time t = t (< 0), that is 



5 These assumptions become worse as the energy of dec rease s below < lOGeV. We discuss the validity of our results by comparing 
with more realistic calculation using the GALPROP code I172H in Sec lIVI 
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The flux at r is given by $(e,r) = (c/4n)f(e, r)[m~ 2 ■ s _1 ■ sr -1 ■ GeV -1 ]. 

Now, in order to estimate the observed flux, we have to specify the cooling function P(e), diffusion coefficient D(e) 
and injected energy spectrum Q(e). First, we formulate the cooling function including the KN effect. Following 
the equation (5) in [92|], the energy loss rate of the e^ 1 including the KN effect is written as 

de P hUtot(e P h)/KN ^ TOe 2 P c 4 ^) 



P(e) 



de 
"dt 



-arc 



r>2 

■°ISM 

8tt 



(53) 



Here or = 6.62 x 10 _25 cm 2 is the Thomson scattering cross section and e p h are the energy of the background photon. 
-Bism is the magnetic field strength in the ISM where we set Bi SM = lfiG. /kn is the KN suppression function which 
is explicitly shown in [l6C)j |. 



Mb) 



63 1 



where 

9(b) = 

and Li2 is the dilogarithm 



b) 
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96- 



(1 + 6)' 



-2 + 2Li 2 (-6) 



Li 2 (z) = 



ln(l - s)ds 



The ISM photons consists of the stellar radiation, reemitted radiation from dust, and CMB, 

Wtot = "star + "dust + "CMB • 



(54) 



(55) 



(56) 



(57) 



Here we model the interstellar radiation field using the results of the GALPROP code [161]. Fig|3] shows the ISM 




energy [eV] 



FIG. 3: The energy density of the ISM photon field at 8kpc from the center of our Galaxy [16 



radiation field energy density e p h X Utot(eph) at ~ 8kpc from the center of our Galaxy. Following the formulation 
above, we numerically calculate the e ± cooling function including the KN effect. Figf4] shows the cooling function 
for e ± with or without the KN effect. The solid line shows the function P(e) in Eq. ([53| . The dotted line shows the 
cooling function when we set /kn = 1- We can see that the KN effect becomes relevant for e > ITeV. 

Second, we formulate the diffusion coefficient. As the diffusion coefficient D(e) for the propagating through our 
Galaxy, we use an empirical law given by the boron-to-carbon ratio observation, that is 

D(e) = D (l + -^-Y . (58) 



3GeV/ 

Here D Q = 5.8 x 10 28 cm 2 s-\ S = 1/3 (H- 

Finally, we assume that the intrinsic energy spectrum at the source is described by the cutoff power law, that is 



Q(e,to,i) = Qq€ "exp 



£cut 



t -t 



(59) 



Here r is the lifetime of the source and to is the time when the leave the source. Then substituting Eq. ([53| . (|58|) 
and (|59p into Eq. (|5ip . we can get the observed electron distribution function f{e,r,t) from a pulsar which is located 
at the distance r from the solar system and born at t = t. 
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10° 



energy [GeV] 



FIG. 4: The cooling function P(e) for e at ~ 8 kpc from the center of our Galaxy. 



B. e distribution function from multiple sources 

Here we consider the e ± distribution function from multiple sources. As we show in Sec lII CI there should be 
multiple pulsars which contribute to the observed e ± flux. 

To calculate the distribution function from multiple sources, we integrate Eq. (|51|) for the pulsar birth time t and 
the pulsar position r, taking into account the birth rate of the pulsars. Then the observed e ± distribution function is 



F(e) 



dt 



»-dif(e.£f) 



2irrdr ■ a ■ ?7wD/(e, r, t). 



(60) 



Again t/wd is the merger rate of the double degenerate WD binary, and a is the fraction of the mergers resulting 
in WD pulsars. We take the lifetime of the WD pulsars t — -t^d as the lower limit of the time integral. We have 
comfirmed that the following results does not depends on this limit as long as it is smaller than — 7\vd- As the upper 
limit of the space integral we take the diffusion length r<iif(e, e t -), which is defined in the same way as Eq. (p)0")) . Through 
the distance rdif(e, e t -), the energy of the propagating e ± changes from to e. 

Since Eq. (|60]) is the mean value, we also estimate the standard deviation of the calculated energy spectrum, that is 



(SF) 2 = 



r r dit {e,e { ) 

2 - 1 di 



2irrdr ■ a ■ r] WD f 2 - Nf. 



2 

avc 



-T W D JO 

where N is the number of the pulsars in our Galaxy that are the source of observing , that is 

N = I di 2-Krdr ■ a ■ ?7wd , 

-TWD J 



(61) 



(62) 



and /avc = F{e)/N is the averaged spectrum per pulsar. We should note that the integral of Eq. (f6"Tj) conta ins a 
serious d ivergence at i = because of the large but improbable contributions from very young and nearby sources jl!5l . 
Il38l4l40| . Here we follow Ptuskin et al (2006) and set the cutoff parameter as 

£ c (e) - -(^wd ■ aD(e))- 1 / 2 , (63) 

which approximately corresponds to the birth time of the newest pulsr that contributes e with energy e. 



Results 



Here we consider two types of models. FigJS] shows the WD pulsar dominant model. In the left panel, the e^- 
flux from multiple WD pulsars (thin solid line) is shown with the standard deviations (thin dashed line), background 
flux (dotted line) and total flux (thick solid line). For each WD pulsar, we set the cutoff energy of the injection 
spectrum e cu t ~ ITeV fEq. (|59p ). intrinsic spectral index v — 1.9, lifetime 7wd ~ 10 9 yr, total energy for each 
~ 10 50 erg, merger rate of double degenerate WD pulsar binaries r\ = 10 _5 yr _1 kpc~ 2 and probability of forming 
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FIG. 5: WD pulsar dominant model : The left panel shows the energy spectrum of the e including the contributions of WD 
pulsars (thick solid line). In this model, we set the cutoff energy of the injection spectrum for each WD pulsar e cut ~ ITeV. 
The average flux (thin solid line) , flux with the standard deviation (thin dashed lines) and background (dotted line) are shown. 
For the background flux, we adopt the fitting function in Baltz & Edsjo (1999) [o^lwith an exponential cutoff for the primary 
electron flux at 5TeV, which is similar to that shown in Aharonian et al. (2008) Tfl ■ We assume that each WD pulsar emits 
the same amount of e . We set the total energy ~ 10 50 erg for each WD pulsar, intrinsic spectral index v — 1.9, lifetime of WD 
pulsars ~ 10 9 yr and birth rate in our Galaxy a ■ ijwd ~ 10~ 7 yr _1 kpc -2 . The right panel shows the positron fraction resulting 
from the average spectrum (solid line) with the dispersion (dushed lines) and background (dotted line), compared with the 
PAMELA data. The background contribution begins to rise around ~ 3TeV since we set the exponential cutoff only for the 
primary electron background, not for the secondary e background. Note also that the solar modulation is important below 
~ lOGeV. 



WD pulsars a = 0.01, which means that the birth rate of WD pulsars in our Galaxy is - 10" 7 yr _1 kpc~ 2 . Tne 
left panel of FigjS] includes the observational data of cosmic ray electrons plus positrons given by the balloon and 
satellite experiments, ATIC /PPT-BETS /Fermi uiHfiL and also the data of ground-based air Cherenkov telescopes, 
H.E.S.S/KASKADE/GRAPES/CASA-MIA [fl I a l 16 2^16^ . For KASKADE/GRAPES/CASA-MIA, the plots show 
the observed flux of the diffuse gamma rays. Since a gamma-ray entering into the air first produces an e^- pair to 
begin a cascade, its shower will look very similar to that of an of equivalent energy [la ]. Thus we presume these 
date as the upper limits on the flux. H.E.S.S electron data are also partly contaminated with photons. Therefore, 
a viable model should not significantly overshoot the points. The background flux consists of the primary electrons 
which is conventionally attributed to the SNRs and the secondary e^- produced by the hadron interaction between 
cosmic ray protons and the interstellar matter, and successive pion decays. For the secondary flux, we adopt the 
fitting function in Baltz & Edsjo (1999) (3^ [^ . [TToj . For the primary electron flux, we also refer Baltz & Edsjo (1999) 
but with an exponential cutoff at 5TeV 6 , which is similar to that shown in Aharonian et al. (2008) 17]. Our result 
fits well the observational data of H.E.S.S and Fermi. 

The right panel of FigJS] shows the positron fraction using the same parameters as the left panel. The results show 
that the observed positron excess can be explained by considering only the contribution from multiple WD pulsars, and 
the positron fraction is expected to drop at around the WD pulsar cutoff energy ~ TeV. The background contribution 
of the positron fraction begins to rise around ~ 3TeV since we set the exponential cutoff only for the primary electron 
background, not for the secondary background. This treatment is appropriate since the abundance of cosmic ray 
protons is observationally robust in this energy range and so is the amount of the secondary background. As we 
discuss in Sec lIVl our calculations become less reliable below < lOGeV since we neglect the anisotropic effects during 
the diffusion in the Galactic disk. Note that in these energy range, the solar modulation is also relevant. 

Figj6] shows the WD and NS pulsar mixed model. In the left panel, the thin solid line shows the flux from 
multiple WD pulsars which have the total energy for each ~ 5 x 10 49 erg, cutoff energy of the injection spectrum 



We also reduce the flux by 30% since the fitting function of Baltz & Edsjo (1999) provide larger flux than the data even without other 
contributions. 
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FIG. 6: WD and NS pulsar mixed model : The left panel shows the energy spectrum of the e (thick solid line) including the 
contributions of WD pulsars (thin solid line) with the same parameters as Fig[5]except for the total energy for each ~ 5 x 10 49 erg 
and cutoff energy of the injection spectrum e cu t ~ lOTeV. The contribution of multiple NS pulsars (dotted dash line) with 
the total energy ~ 10 48 erg for each NS, cutoff energy of the injection e cut ~ ITeV, lifetime ~ 10 5 yr, birth rate in our Galaxy 
~ lCP 5 yr _1 kpc -2 are also included. For both fluxes of pulsars, the standard deviations are shown. The total e flux and its 
deviation are the thick solid and dushed line, respectively. The right panel shows the positron fraction. 



£cut ~ lOTeV and the same value for other parameters as Figj5j The difference of the e cu t means, in our WD pulsar 
model, the difference of the multiplicity Ai, the magnetic field strength B p , the angular freauency f2 and the radius 
R according to Eq.(|20p. The dotted-dash line shows the e flux from multiple NS pulsars with the total energy 
~ 10 48 erg for each, cutoff energy of the injection energy spectrum e cu t ~ ITeV, lifetime ~ 10 5 yr, birth rate in our 
Galaxy ~ 10 — 5 yr _1 kpc _ . The standard deviation of the e ± energy flux from the WD pulsars is relatively small 
compared with the NS pulsars. This is because the larger abundance of WD puslar are expected as we discussed 
in the previous sections. The dotted line shows the same background contribution as FigJSJ The total flux and the 
deviations are shown by thick solid line and thick dushed line, respectively. It is shown that the excess in the e ± flux 
in the range lOOGeV < e < ITeV is explained by the multiple NS pulsars. By considering the contribution of multiple 
WD pulsars, the smooth "double bump" are forme d in the e nergy spec trum around ITeV and lOTeV, which can be 
observable by the future experiments like CALET fulfill] and CTA [l20l |. 

The right panel of Fig|5] shows the positron fraction for the mixed model. The observed positron excess can be 
explained and, in this case, there will be no flux drop around ~ TeV in contrast to Fig. [5] 



D. Observed WD pulsar candidates 



Finally, in this subsection we give two interesting examples of the observed WD pulsar candidates, AE Aquarii 
and EUVE J0317-855. So far, a few thousan d WDs h ave been discovered, and the magnetic field and the rotational 
period have been detected for some of them 165- 168]. Forthcoming experiments like ASTRO-H 169] will find more 



magnetized and rapidly spinning WDs, which will reveal the detailed characteristics of such WDs. Then, we will 
know whether sufficient amount of WD pulsars exist in our Galaxy or not. 

AE Aquarii 

AE Aquarii is a magnetized cataclysmic variable, that is, consisting of a primary WD and a spectral type K5V 
main sequence star, located at ~ lOOpc from the solar system. The primary WD has spin period ~ 33 s, which i s 
identified by the approximately sinusoidal profiles of the observed emissions at energies below ~ 4keV jl52l . Il53| ] . 
Recently Suzaku satellite discovered that AE Aquarii shows hard X-ray sharp pulsations at t he p eriod consistent 
with its rotation [1361 ] . Also TeV gamma emissions during the optical flares were reported [154 Il55l ] , although there 
have been no detection since then. The primary WD is spinning down at a rate ~6x 10~ 14 sec sec -1 , implying the 
spin down luminosity ~ 10 33 erg/sec, which is three orders of magnitude larger tha n the UV to X-ray emissions. The 
magnetic field strength inferred from the spin down luminosity is ~ 5 x 10 7 G PH. 

Since the AE Aquarii is an accreting binary system, the density of the plasma surrounding the primary WD may 
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be much higher than the GJ density. However, both theoretical (157j and observational works suggest that the rapid 
rotation and strong magnetic field produce a low-density region around the WD, and the particle acceleration by 
the same mechanism as spin-powered pulsars could be possible. The parameters of AE Aquarii satisfy the condition 
Eq.JUJl, above the death line of WD pulsars (FigHJ). 

EUVE J0317-855 (RE J0317-853) 

EUVE J0317-855 is a hydrogen-rich magnetized WD discovered by ROSAT and EUVE survey [HI, Il59j |. By 
analyzing the photometric, spectroscopic and polarimetric variations, EUVE J0317-855 is shown to rotate at the 
period ~ 725s, which is one of the fastest isolated WDs, and the dipole magnetic field is ~ 4.5 x 10 8 G. EUVE 
J0317-855 have a DA WD companion which is located at > 10 3 AU from EUVE J0317-855. Because of the large 
separation, there sup pose to be no interaction between the two WDs. By analysing the emission from the companion, 
Barstow et al (1995) |l58| noted that EUVE J0317-855 is located at ~ 35pc from the solar system, and the mass is 
1.31-1.37M© which is relatively large compared with the typical WD mass ~ 0.6M Q . Its rapid ro tatio n and large mass 
suggest that EUVE J0317-855 may be the outcome of a double degenerate WD binary merger |l59| . Relevant pulse 
emission from EUVE J0317-855 has not been observed yet, which may suggest that the creation and acceleration 
does not occur. When we put the parameters of EUVE J0317-855 on Fig[TJ it comes below the death line, which is 
also consistent with the observation. 



IV. SUMMARY AND DISCUSSION 

We have investigated the possibility that WD pulsars become a new TeV source. We have supposed that a fair 
fraction of double degenerate WD binaries merge to become WD pulsars, and these WDs have the magnetospheres 
and pulsar wind nebulae. The pair creation in the magnetospheres and their acceleration and cooling in the wind 
nebulae have been discussed, and we have found the following. 

1. If a double degenerate WD binary merges into a maximally spinning WD, its rotational energy will become 
~ 10 50 erg, which is comparable to that of a NS pulsar. Also the birth rate ~ 10~ 2 -lCU 3 /yr/galaxy is similar to 
the NS case, which provides the right energy budget for cosmic ray e ± . 

2. Applying the theory of NS magnetospheres, we give the e ± pair creation condition ("the death line") for WD 
pulsars. Since our fiducial parameters of WD pulsars meet the condition, the WD pulsars are eligible for the 
e ± factories. The death line is consistent with the observations for some WD pulsar candidates. 

3. By assuming the energy equipartition between e^- and magnetic field in the wind region, we have shown that the 

produced in the WD pulsar magnetosphere can accelerate up to ~ lOTeV when the WD pulsar has a rapid 
rotation (P ~ 50s) and strong magnetic fields (B ~ 10 8 G) and the e ± multiplicity is not so large (M ~ 1). 

4. In contrast to the NS case, the adiabatic energy losses of e ± in the pulsar wind nebula region are negligible 
in the case of the WD pulsars since they continue to inject the e ± after the nebula stop expanding. Also the 
radiative cooling of is not so large, and the high energy can escape from the nebula without losing much 
energy. As a consequence, it is enough that a fraction ~ 1% of WDs are magnetized as observed in order for 
the WD pulsars to become the relevant TeV sources. 

Based on the WD pulsar model above, we have calculated the observed flux from multiple WD pulsars in our 
Galaxy. We have solved the diffusion equation including the KN effect, and found the following. 

5. We have shown the two model fluxes. In one model (WD pulsar dominant model), only considering the 
contribution from the multiple WD pulsars, we can explain the reported excess of the flux around lOOGeV < 
e < ITeV and also the PAMELA positron excess. In the other model (WD and NS pulsar mixed model), 
the combination of the multiple WD and NS pulsars can also explain the existent observations, and form the 
double bu mp in the energ y spectrum of e^, which can be a signature for the future observation like CALET 
[HI EH and CTA pof . Since the lifetime of WD pulsars are relatively large, the number of nearby active 
sources can be huge, which give a small Poisson fluctuation in the e^- flux compared with NS pulsars. 

As we have shown, WD pulsars could dominate the quickly cooling above TeV energy as a second spectral bump 
or even surpass the NS pulsars in the observing energy ra nge ~ 100 GeV, pro viding a b ackgr ound for the dark matter 
signals and a nice target fo the future AMS-02 [ll^ . 1ll7j . CALET [HI EDI and CTA [l20|. As the future works we 
should consider other observational signatures than e^- for the coming multi-messenger astronomy era. For example 
we have to consider the radio to 7 ray emission from WD pulsars based on our model. The number of observed 
pulsars in the Galactic disk should be proportional to ~ (number density) x (radio luminosity). Since about ~ 10 3 
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NS pulsars have been discovered by radio telescopes, assuming that WD pulsars can convert the spin down luminosity 
to the radio emission with the same radio efficiency as NS pulsars, the number of WD pulsars which should have been 
already detected by radio observation can be estimated as 

10 3 ( *-Vwd \ f twd\ / Lwd\ ^ 1Q ,_a_\ (g4) 



Vns J \tns J \LnsJ vO.Ol 

Thus O(10) WD pulsars may well be observed as radio pulsar with relatively long period P ~ 50sec. However, since 
the efficiency of the radio emission depends on the detailed situation in the polar cap regions, whether WD pulsars 
have the same efficiency as NS pulsars is highly uncertain at this stage. Other than the electromagnetic emissions, 
double degenerate WD mergers, which w e co nsider as the origin of WD pulsars, is a promising source of the future 



gravitational wave observation by LISA [131| . It is very interesting if we get a strong constraint on the event rate 
of the mergers in our Galaxy by observing the high energy e . In this paper, we consider only merged WDs as a 
source of high energy e emissions. In the single-degenerate binaries, the accretion could induce the rapid rotation of 
the WDs. These accreting binary systems could also become WD pulsars if they have strong magnetic fields as AE 
Aquarii. 




Energy e e [GeV] 



FIG. 7: The comparison of our result and the GALPROP code for e total flux from multiple NS pulsars. 

At the current moment, our model have several crucial assumptions which should be considered more carefully. 
Last of all, we discuss these things and what we should consider in the future works. 

A. We have discussed the death line of WD pulsars based on the simplest polar cap model, considering only the 
curvature radiation for the pair creation photons and the vacuum polar cap gap in which p = 0. It has been 
shown that the inverse Compton scattered photon is important for the pair creation in the pola r cap, and 
the observed death line of NS pulsars is well explained also by the space charge limited flow model jl70l . Il7ll | . 
Especially for the WD, even in the case of f2B < 0, the charged limited flows may exist since the binding energy 
of the ions could be smaller than the thermal energy at the surface. Henc e we have to investigate the death line 
for WD pulsars based on, for example, the Harding & Muslimov model [l70l Il7l| . Also the e ± multiplicity in 
the magnetosphere is crucial for the maximum energy, and we have to calculate it consistently with the polar 
cap model. 

B. There are uncertainties about the accelerating and cooling processes of in the pulsar wind ne bula . Th e 
energy flux of the magnetic field may not be conserved in the wind region like the Crab nebula [l50l Il5lj . 
This time we have assumed Eq. (|19j) for simplicity. Moreover, we have to evaluate more precisely the inverse 
Compton scattering in the pulsar wind nebula as a radiative cooling. In this paper we roughly estimate it to be 
comparable to the synchrotron radiation. We also have to worry about whether the wind mainly consist of the 
e^, which is still under debate even in the case of NS pulsars. 

C. When calculating the flux from the multiple sources, we assume that the source distribution and the 
diffusion process are isotropic. However compact objects like WDs and NSs may distribute more densely near 
the center of our Galaxy. (Also, the large kick which could be given at their birth may affect the spatial 
distribution of the pulsars.) Since the arrival anisotropy can be useful to discriminate the origin of the observed 
e^, we should take into account these anisotropic effects. For the with relatively low energy < lOGeV, the 
inverse Compton energy losses become less important and consequently the propagation range of e^- increases, i.e. 
the anisotropic effects during the propagation, for example the effect of the Galactic disk structure, become more 



18 



prominent. In these low energy range, the public GALPROP code [1721 ] can provide a more reliable calculation 
of the propagation from distant sources arbitrar ily d istributed in our Galaxy. FigJT] shows the comparison of 
our result and the GALPROP code (WEBRUN |l72| ) for the primary e ± total flux from multiple NS pulsars 
with the same parameters as FigO We have confirmed that our result is consistent with the more realistic 
calculation in the high energy region and begin to deviate below < lOGeV. The bump around 0.5GeV in the 
result of the GALPROP code is formed mainly due to the diffusive reacceleration of e ± during the propagation 
in our Galaxy. Note that in this region, the solar modulation is relevant and the uncertainty becomes large. 
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